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ABSTRACT 

We analyzed XMM-Newton observations of two center-filled supernova rem- 
nants, G27.8+0.6 and G28.8+1.5, to search for pulsars/pulsar wind nebulae 
(PWNe) and other types of neutron stars associated with these remnants. We 
discovered a PWN candidate within the extent of the centrally-peaked radio 
emission of G27.8+0.6. The X-ray morphology of the PWN candidate and its 
position with respect to the host supernova remnant suggest that the alleged pul- 
sar might be moving away from the supernova remnant's center with a transverse 
velocity of around 100-200 km s _1 . The majority of the detected X-ray point 
sources in both fields are classified as main-sequence stars, based on their bright 
optical counterparts and relatively soft X-ray spectra. The remaining medium- 
hard and hard sources are most probably either AGNs or cataclysmic variables, 
although we cannot completely rule out the possibility that some of these sources 
are neutron stars. 

Subject headings: SNR: individual (G27.8+0.6, G28.8+1.5) - - SNR: ISM - 
stars: neutron — X-rays: stars 



1. Introduction 

Stellar evolution models predict that the core-collapse supernova explosions produce 
both a supernova remnant (SNR) and either a neutron star (NS) or a black hole (BH). The 
exact fraction of these two types of compact objects is not well constrained, but considering 



1 Based on observations obtained with XMM-Newton, an ESA science mission with instruments and 
contributions directly funded by ESA Member States and NASA. 
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that the majority of SN explosions are believed to be the core-collapse events, a significant 
fraction of SNRs are expected to harbour NSs. A significant number of these NSs are rapidly 
spinning, and are often detected as radio, X-ray or 7-ray pulsars. 

A young energetic pulsar loses most of its rotational energy in form of a pulsar wind 
( Rees fc Gunn|[l~974 ) , which interacts with the surrounding medium and forms a pulsar wind 



nebula (PWN) that is observed as radio and/or X-ray diffuse emission (Gaensler & Slane 



2006). Due to the complexity of this interaction and anisotropy of the pulsar wind, PWNe 



display various shapes including toroidal structures, jets, and sometimes comet ary- shaped 
tails formed when the pulsar is moving with a supersonic speed through the ambient medium 
(Kargaltsev & Pavlov [2008 ) . Hence, to obtain a better census of the NS-SNR associations, 
high-resolution X-ray (and radio) observations are very useful, as they can identify PWNe. 

Recent sensitive, high-resolution X- and 7-ray observations have also revealed unex- 
pected diversity in the populations of NSs associated with SNRs, which, in addition to 
classical radio pulsars, also include anomalous X-ray pulsars (AXPs), soft gamma-ray re- 
peaters (SGRs), compact central objects (CCOs), and perhaps other classes of radio-quiet 
NSs QKaspi fc Helfandl|2002| ) . A class of radio-quiet NSs was predicted by |Shklovskii| ( |1980[ ) , 
who suggested that due to their strong magnetic fields, these NSs will decelerate rapidly and 
become inactive in several hundreds of years after the SN explosion. Alternatively, NSs can 



be born with low magnetic fields ( "antimagnetars" ; e.g., see Halpern & Gotthelf 2010, and 



references therein). The CCO in CasA, revealed in the first-light Chandra image (see also 
Pavlov fc Lu"rla|[20"09l ), is an example of such a source. 



Out of almost 50 NSs possibly associated with SNRs (listed in Table 1 of Kaspi et al. 
20021, -30% were detected as radio pulsars, ~30% as various radio-quiet compact objects 



(AXPs, SGRs, CCOs), while the remaining NS-SNR associations were proposed because a 
diffuse, non-thermal X-ray nebula (PWN), and/or radio synchrotron nebula (plerion) were 



detected within the SNR, suggesting the presence of a pulsar. Kargaltsev &; Pavlov (2008) 



list a total of 18 SNR-PWN-pulsar associations and 14 PWNe in which no pulsar has been 



detected to date. Kaplan et al. (2004, 2006) searched for NSs in a volume-limited (closer 



than 5kpc) sample of radio-shell SNRs, but found no NSs down to the luminosity limit of 
typical NSs of the same ages as the observed SNRs. To search for compact stellar remnants 
and/or PWNe in all center-filled and composite SNRs within 5 kpc, a snapshot survey has 
been carried out with XMM-Newton. We present the analysis of two such SNRs, G27.8+0.6 
and G28.8+1.5. 

The center-filled SNR G27.8+0.6, with a relatively low surface brightness, was detected 
in single-dish pointed observations at 2.69, 4.75 and 10.2 GHz (Reich et al. 1984). The 



central part is strongly polarized, indicating a non-thermal origin and supporting the SNR 
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classification. Reich et al. (1984) estimated the age and distance of the SNR to be 35-55 



kyrs and 2-3 kpc, respectively, adopting the evolutionary model for center-filled SNRs of 



Weiler k Panagia (1980). 



G28.8+1.5 was detected by Schwentker (1994) in ROSAT observations as a diffuse X-ray 



emission covering a region of 80' x 50', and identified by non-thermal radio emission as a 
~30-kyr-old SNR at a distance of ~4 kpc. Schwentker (1994) also reported a 5.45 s period 



pulsations from a point source (which is not in the XMM-Newton field of view) in the vicinity 



of this SNR. However, this result has not been confirmed later (Song et al. 2000) 



In a search for compact objects in these two center-filled SNRs, we analyzed all X- 
ray sources detected by XMM-Newton in their fields. To classify the detected objects and 
select PWN/NS candidates, we used the inferred X-ray properties and cross-correlation with 
optical, infrared, and radio data. In Section|2]we summarize our data analysis and present the 
compiled catalogs of the detected X-ray sources in both fields (Section |2TT), and also a search 



for their counterparts at other wavelengths and the proposed classifications (Section 2.2). 
Finally, we discuss PWN/NS candidates in Section |3| 



2. Observations and results 

G27.8+0.6 and G28.8+1.5 were observed in the Full Frame mode using medium filters. 
Table [T] lists the observation IDs and dates, the total exposure times and the exposure times 
after high-background screening. The data analysis was performed using the XMM-Newton 
Science Analysis System (SAS) v7.1.0. We produced PN, MOS1, and MOS2 event files using 
the SAS tasks epchain and emchain, and then excluded the times with high background 
using a threshold count-rate of 1.5 cts s~ x for all three detectors. 

To search for X-ray sources in the two SNR fields, we produced PN, MOS1 and MOS2 
images in five energy bands: 0.2-0.5 keV, 0.5-1.0 keV, 1.0-2.0 keV, 2.0-4.5 keV, and 4.5-12.0 
keV (bands Bl to B5, respectively). This is a standard band selection for the source detec- 
tion procedures in SAS, for which the source significance (detection likelihood) calculation 
is well established. For the PN images in band Bl we selected only single-pixel events (PAT- 
TERNS), while single and double-pixel events (PATTERN<4) were included in the other 
energy bands. Such a PATTERN selection was motivated by significantly higher noise from 
particle background in the softest PN band. We selected 1-4 pixel events (PATTERN<12) 
in all five bands of the MOS images. For each image we produced the exposure map corrected 
for vignetting, the background and sensitivity maps, and the corresponding detector mask 
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Fig. 1 — Left: VLA images of G27.8+0.6 at 1.4 GHz (top) and G28.8+1.5 at 0.325 GHz 
(bottom). Both images are extracted from MAGPIS (White et al. 2005). The radio beams 
(~ 6" x 5" for G27.8+0.6 and ~ 66" x 45" for G28.8+1.5) are shown in the lower left corner of 
each image. Right: The broad band (0.2-12.0 keV) gray scale XMM-Newton EPIC images 
of the observations 0301880401 (top) and 0301880901 (bottom), combining PN, MOS1, and 
MOS2 data. The images are spatially aligned with the corresponding radio images, and 
smoothed with a Gaussian of FWHM 5". The intensity scale is linear from to 2.6943 ct s -1 
pixel -1 . The pixel size is 4". The circles (radius 36") represent the sources detected above 
the likelihood threshold of 7. The numbers correspond to the SRC IDs from Tables [3] to [6j 
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to define source detection areas. These maps were used as inputs for our source detection 
procedures in each SNR field. 



2.1. The X-ray source catalogs 

To detect the X-ray sources in the observed fields, we run the standard SAS source 
detection procedures searching simultaneously in a total of 15 images for each field (the 5 
energy bands Bl to B5 in each of the 3 EPIC cameras). The SAS task emldetect was used 
to determine the parameters of the detected sources (significance, count-rates, fluxes, and 
hardness ratios) by using a maximum likelihood fit. We selected the likelihood threshold 
value of mlmin=7, which allows detection of even very weak sources (S/N^3— 4; see the SAS 
manual) . 

To convert the vignetting-corrected count rate to flux for each source in each energy 
band, we estimated the energy conversion factors for the PN and MOS medium filters using 
W3PIMMS0 (Table The cataloged values are calculated for a power-law model with a 
photon index of 1.7 absorbed by the A^h columns of 1.5 and 2 x 10 22 cm -2 for G27 and G28, 
respective!}]^} This model fits well the majority of the detected sources, i.e, intrinsically 
medium-hard sources, although we found that the inferred fluxes are not very sensitive to 
the spectral model selection. For example, for soft sources (e.g., best-fitted by a black-body 
model) our cataloged flux values are overestimated by only ~10%. On the other hand, 
for significantly harder sources (for example, a power- law with a photon index of ~1) the 
cataloged flux values are underestimated by approximately 30%. The calculated fluxes are 
actually much more sensitive to the assumed absorption columns. For example, nearby stars 
would be much less absorbed, but even if we assume an A% column of only 5 x 10 20 cm -2 , 
the flux values listed in Tables [3] and [5] would be overestimated by approximately a factor of 



2, which, however, does not affect the source classification significantly (Section 2.2) 



After compiling the list of X-ray sources detected above the 3a — 4a threshold in both 
SNR fields, we compared the X-ray positions with optical and near-infrared (NIR) data 
to determine possible offsets of the two XMM-Newton pointings. We cross-correlated the 



X-ray positions with optical (USNO-B1; Monet et al. 2003) and NIR (2MASS; Skrutskie 



et al. 2006) catalogs using a distance of < 3x the combined X-ray and optical/infrared 



positional error to associate the optical/infrared and the X-ray source. We used the identified 



2 ht tp : / /heasarc . gsfc .nasa. gov / Tools / w3pimms . html 

3 These absorption columns are equal to the Galactic Hi columns measured in the directions of the two 



SNRs (Table [2]). 



-6- 



optical/infrared counterpart candidates to calculate mean offsets in RA and DEC, but did 
not notice any systematic shift (within the uncertainties) with respect to the optical/infrared 



data. A standard la (residual) systematic positional error of 1'.'5 (Watson et al. 2003) was 



then added quadratically to the statistical positional uncertainties of the cataloged sources. 

We detected 11 point X-ray sources in the field of G27.8+0.6, and 28 sources in the 
G28.8+15 field. The results are summarized in Tables [3] to [6] Tables [3] and [5] list the X-ray 
properties of these sources. The source identification number is listed in column 1 (the same 
source number is also used in Figure [T]), while the position and the total positional error, 
which includes the statistical and systematic errors, are listed in clumns 2 and 3, respectively. 
The source detection likelihoods in all 15 images are added and transformed to the equivalent 
single-band detection likelihood, which is listed in column 4. Columns 5 and 6 list the total 
number of the accumulated counts and corresponding error, and the total flux and error in 
the 0.2-12.0 keV band. The fraction of the soft (0.2-4.5 keV) flux, used to calculate the 
ratio of the X-ray to optical flux (see Section 2.2[ ), is listed in column 7. The hardness ratio, 



defined as HR = (-R 2 -i2keV - -Ro.2-2kcv)/(-R2-i2keV + -Ro.2-2kev) (where R is the count rate) 
and errors, all from the combined PN, MOS1, and MOS2 data, are listed in column 8. Visual 
inspection of the broad-band images revealed a marginaly extended source (src. #1) in the 
field of G27.8+0.6. We discuss this source in Section [3] 

Tables [4] and [6] list the results of the cross- correlation with the USNO-B1, 2MASS, and 
radio (NVSS, [Condon et al.||1998~ ) catalogs, which are described in detail in Section 2.2 The 



detected USNO-B1/2MASS counterparts are listed in column 2, together with the offsets 
from the X-ray position (3), and the corresponding X-ray to optical flux ratio (4). The 



tentative source classification is given in column 5, and we discuss it further in Section 2.2 
while column 6 lists the probability of the optical/NIR association. 

The point-source sensitivity limit of the XMM-Newton observation in the field of G27.8+0.6 
is about 1.5 x 10~ 14 erg cm -2 s _1 [corresponding to an unabsorbed luminosity of (1-3) xlO 31 
erg s _1 at the estimated SNR distance of 2-3 kpc and Nn = 1.5 x 10 22 cm -2 ], and slightly 
deeper across G28.8+1.5 [(1 — 1.5) x 10~ 14 ergs cm -2 s _1 or (2 — 6) x 10 31 erg s _1 at the 
distance of 3-5 kpc and = 2 x 10 22 cm" 2 estimated for this SNR]. The brightest source 
in the G27.8+0.6 field (~ 2 x 10~ 12 ergs cm -2 s -1 ) is an order of magnitude brighter than 
any of the other detected sources, while in the second field the measured absorbed fluxes 
range from 1 x 10 -14 to 3 x 10 -13 erg cm -2 s _1 . 

Figure [l] shows the radio images of the two SNRs, and the broad band (0.2-12.0 keV) 
images of both fields, combining PN, MOS1, and MOS2 data. The detected X-ray sources 
are shown as circles in each field. The top left panel shows the VLA (1.4 GHz) image of 
G27.8+0.6. Its bright central region is clearly visible as well as a larger (30' x 50') faint 



-7- 



extended component. The fainter extended component is also seen at lower frequencies in 



the single-dish radio observations (Reich et al. 1984). The bottom left panel shows the 



region of G28.8+1.5 observed by VLA at 0.325 GHz. Both VLA images are extracted from 



MAGPIS (White et al. 2005). Sources 3 (in G27.8+0.6) and 6 (in G28.8+1.5) appear to be 



well centered within the bright radio emission regions (see also Section [3]) . 

Since the XMM-Newton observations of the G27.8+0.6 and G28.8+1.5 are very short, 
the small number of accumulated counts limits spectral analysis, even for the brightest 
sources. However, the hardness ratios suggest that we detected different source populations: 
the sources exhibiting intrinsically hard spectra (e.g HR> 0), and several extremely soft 
sources with most of the X-ray emission below 0.5 keV (e.g HR~ —1). 



2.2. Cross-correlation and classification of X-ray sources 

We cross-correlated all the X-ray sources listed in Tables [3] and [5] with the USNO-B1, 
2MASS, and NVSS catalogs. 

We found USNO-B1 counterpart candidates for each of the 28 X-ray sources detected in 
the G28.0+1.5 field, and all except two sources (src. #1 and #7) in the field of G27.8+0.6. 
Although we find most of the candidate optical counterparts within the ~ 2a positional 
uncertainty of the X-ray sourc^j there is a high probability of a chance coincidence, due 
to a broad PSFs of the XMM-Newton detectors and a high density of stars in the Galactic 
plane. 

Based on the optical source density (p = 0.0080 and p = 0.0077 field sources in a 1" x 1" 
box in the field of G27.8+0.6 and G27.8+1.5, respectively), we calculate the probababilit}|^] 
that each listed optical source is the true counterpart (Tables [1] and [fT 



To classify the sources we also used the X-ray-to-optical flux ratio calculated as log(/ x / /< 



opty 



log(/ x ) + (m-B2 + "^R2)/(2 x 2.5) + 5.37, following Maccacaro et al. (1988). In this formula 



f x is measured in the first four bands (0.2-4.5 keV) because the source classification scheme 



of Maccacaro et al. (1988) is based on the relatively soft X-ray band (0.3-3.5 keV) of the 
Einstein Observatory. We used the average of the B2 and R2 USNO magnitudes instead 
of (mostly unavailable) V magnitudes, but the difference between these magnitudes for the 
main sequence stars should be less than ±0.1 (?). We have verified this for several stars in 



4 All the optical candidate counterparts are within the 3a of the X-ray positional error. 

5 We estimate it as the probability of finding zero field sources in the circle of radius d, P — exp(— pud 2 ), 



where d is the distance between the optical and X-ray positions. 
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Fig. 2. — The combined MOS1+MOS2 image (left) and DSS image (right) showing the 
region around the brightest point source (src. 1) in the field of G28.8+1.5. The black 
circle in the DSS image shows the position of the X-ray source as determined by our source 
detection procedure (R.A = 18:39:03.041, Decl. = -02:42:01.94). The radius of the circle is 
4.5" (3a positional error). The two USNO-B1 sources, which are most probably the optical 
counterparts of the detected X-ray emission, are also marked in the DSS image (see text for 
more detail). The 3a error circle of the cataloged X-ray position is also indicated on the 
MOS image (as the white circle). 
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our catalogs, for which the USNO and the position of the optical source with the cataloged 
V magnitude are sufficiently close (i.e., the three magnitudes are most likely measured from 
the same source). However, the difference between the averaged USNO and V magnitudes 
for sources classified as AGN/CVs seems to be slightly larger (up to 0.5) but still not large 
enough to change our source classifications. 

The cross-correlation with the 2MASS catalog yielded 9 and 23 candidate counterparts 
for the G27.8+0.6 and G28.8+f .5 field, respectively. We found only one radio counterpart - 
the 3.6-mJy NVSS source 183932-023919 - coinciding with the X-ray source #27 from the 
G28.8+1.5 field. 

We used the multi-wavelength data to classify the cataloged sources. In particular, we 
used the X-ray hardness ratios and optical and NIR properties to identify main-sequence 
stars and other source classes in our sample. Main-sequence stars are expected to have 
significantly smaller values of log(/ x / / op t)- Depending on the stellar type, the log(/ x / / op t) for 
stars is in the range -5 to -0.5, while galaxies, AGNs, X-ray binaries (XRB) and cataclysmic 
variables (CVs) are expected to have significantly stronger X-ray emission comparing to 



optical (e.g., —1 < log(/ x // opt ) < 1 for AGNs; Maccacaro et al. 1988). We also used NIR 



colors, in addition to the optical data, to confirm the stellar classifications (e.g., see Table 



6 in Finlator et al. 2000). We used these selection criteria and tentatively classified the 



majority of the detected sources in both SNR fields. 



2.2.1. Sources m the G27.8+0.6 field 

In the field of SNR G27.8+0.6 we classified 8 sources as main sequence stars (Table 
|2J) based on their log(/ x // op t) and relatively soft X-ray spectra. The NIR colors support 
the proposed spectral classification, suggesting that the majority of these objects are nearby 
stars. 

There is an optical counterpart candidate within the 3a error circle of source #3. Al- 
though the X-ray-to-optical ratio suggests a stellar classification, the distance between the 
optical and X-ray position is almost 6", making the chance coincidence very likely (see Ta- 
ble |1J). In addition, the X-ray spectrum of this source is significantly harder than the typical 
stellar spectra. We identified a 2MASS source closer to the X-ray position (off by ~2"), 
which suggests that the source is either an AGN, CV or a faint X-ray binary. 

We note that the probability of a chance coincidence between the source #8 and its 
possible optical counterpart is also very high (about 40%). However, this source, which 
exhibits a soft X-ray spectrum, is detected at a large offset from the center of the field of 
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view, which reduces the accuracy of the X-ray position. 

Two sources (src. #1 and #7), which display relatively hard X-ray spectra in compar- 
ison with identified field stars, do not have any cataloged optical or infrared counterparts. 
While one of these objects (src. #7) is detected as a weak X-ray source (S/N ~ 4, see 
Table [3]), source $T is the brightest source in the field. Furthermore, the morphology and 
spectral properties of this source suggest that this object could be a NS created in the ex- 
plosion of G27.8+0.6 and moving away from its birth site. We discuss this source in more 
detail in Section [U 



2.2.2. Sources in the G28. 8+ 1.5 field 



All of the 28 sources detected in the G28.8+1.5 field have optical candidate counterparts 
within the 3cr X-ray error circle. Based on the X-ray, optical and NIR data, we classified 
11 sources as stars, while 3 objects are stars or galaxies, according to their NIR colors. The 
remaining 14 sources have significantly larger X-ray-to-optical flux ratios and harder X-ray 
spectra. The surface number density of these hard sources is approximately the same as 



that of the background extragalactic sources (Ebisawa et al. 2005), suggesting that most of 



them are AGNs, although some of them might be either faint X-ray binaries, or, more likely, 
cataclysmic variables, which are considered to be prime candidates for faint Galactic hard 



X-ray sources (e.g., Ebisawa et al. 2005, and references therein). 



With approximately 400 accumulated counts in the PN data, source #1 (Table |5J) is 
the only source in the field of G28.8+1.5 sufficiently bright for spectral analysis. We found 
the best fit model for this source to be an absorbed MEKAL (kT = 0.7 — 0.9 keV) plus 
a hard component represented by a power-law model (r = 2.2 — 2.3). The low value of 
the absorption column (below 8 x 10 20 cm -2 ) indicates that the source is much closer than 
the SNR. We found that the flux obtained from the spectral fit, fx = (0.6 — 1.3) x 10~ 14 
ergs cm -2 s" 1 (0.2—12 keV), is slightly lower than the cataloged value (which assumes a 
PL model with a significantly larger absorption column). We also obtained an acceptable 
fit for this source using a two-temperature plasma model (kT ~ 0.1 and kT ~ 0.8 keV), 



typical of active late-type stars (e.g., Giidel 2004), but could not find an acceptable fit for 
a single-component or any combination of the power-law or blackbody models. However, as 
Figure [^indicates, there are probably two X-ray sources (most likely stars) at this position, 
which are difficult to resolve with XMM-Newton. 
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3. PWN/NS candidates 

3.1. A PWN candidate in G27.8+0.6 

No extended diffuse X-ray emission coming from the SNR shell or its central regions is 
visible in the EPIC images of G27.8+0.6. This is not surprising because after filtering the 
data affected by flaring the remaining good-time exposure was only about 4.9 ks for PN, 10 
ks for MOS1, and 8.4 ks for MOS2 (Table [I), making it very difficult to detect any faint, 
large-scale emission above the background. 

However, we have detected a bright, possibly extended source offset by ~20' (or ~12 
pc at the SNR's distance of 2-3 kpc) from the presumed SNR's center (Figure [3]). Although 
the source was not flagged as extended in our source detection procedures, the smoothed 
(Figure [3]) and binned (Figure [1]) images suggest that there is a small elongation in the North 
direction and perhaps some weak emission to the East from the point source. Although the 
number of accumulated counts is small (see Table [3j because of the short exposures, the 
extended structures are clearly visible, not only in the combined data but also in individual 
PN and MOS images. We note that there is a CCD gap and a bad column across the source 
in the PN data, which cuts off a significant fraction of the point source, and perhaps also a 
part of the candidate nebulaj^j Since the source is detected 10' off-axis, the EPIC PSFs are 
elongated in the North-South direction. However, the entire source extent is not likely to be 
due to the PSF distortion^} In addition to possible extended emission in the vicinity of the 
point source, there are several relatively bright blobs visible in the Southwest and Northwest 
directions, forming a triangular-like shape with the point source at the vertex. These blobs 
(with a signal to noise of ~3-4, estimated in an aperture of 40" x 40" in the 0.7-10 keV, 
measuring the background from the nearby area) might be part of extended emission at a 
distance of ~ 5' — 10' from the source. Figure [4] (right) shows a NIR (2MASS) image of 
the same region with no optical counterpart at the position of the X-ray source. We have 
searched optical catalogs and found no sources within 15" radius with / < 20 mag (and down 



to the limiting magnitude of the USNO-B1 catalog, M v ~ 21; Monet et al. 2003). 



To determine the spectral properties of our PSR candidate, we extracted the source 
spectra from the PN and MOS data and group them with a minimum of 10 counts per bin. 
However, the small number of the accumulated photons (~60 in PN and slightly more in 



6 We have excluded all border pixels from these images (by selecting FLAG=0), to make sure that no 
artificial structures are produced in the vicinity of the source. 

7 For example, compare with the PSF of source 2 in the G28.6+L5 field, which, at approximately the 
same off-axis angle, is less elongated. 
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M0S1+M0S2 after subtracting the background) in our short-exposure observations, and a 
large (~50%) background contribution, did not allow us to constrain well the fitting param- 
eters. 

We have obtained an acceptable fit (x 2 /d.o.f = 0.84/16) for an absorbed power-law 
model with the photon index T = 1.2 ± 0.8 and the hydrogen column density Ah = (3 ± 2) x 
10 22 cm -2 , consistent with the Hi value from Dickey & Lockman (1990). Using this model, 



we estimated the source flux to be aboutn 4 x 10 13 erg s 1 cm 2 in the 1.5-10 keV band, 



implying a luminosity range of Lx 
kpc. 



(2 — 5) x 10 32 erg s 1 for an assumed distance of 2-3 



The right panel of Figure [3] shows the radio map of G27.8+0.6 and the position of 
the PWN candidate with respect to the SNR center. The two arrows (also shown in the 
EPIC image on the left) indicate two possible velocity directions of the pulsar candidate, 
based on the PWN shape. If the blobs extending westwards are part of a "tail" formed by 
a bow-shock, due to the pulsar's supersonic velocity, then, according to the existing MHD 
models (e.g., Bucciantini et al!]|2005 ), the pulsar travels in the direction of the arrow marked 
2, i.e., away from the center of the brightest part of the extended radio emission, which 
probably coincides with the site of the SN explosion. The extended emission to the North 
could be then explained as equatorial emission (a torus; Weisskopf et al.pOOO ). However, we 
cannot exclude the possibility that the pulsar candidate is moving in the direction of arrow 
1 (i.e., the "tail" extends to the North, and the blobs are not associated with the PWN 
candidate). Although, in this case, the pulsar's velocity would be nearly perpendicular to 
the expected radial direction, similar behavior was also observed in other SNR-PWN systems 
and explained by a high initial velocity of the NS progenitor, which caused the SN explosion 
offset from the center of the wind-blown cavity and the remnant's geometric center (e.g., 



similar to the PWN in the IC443; Gaensler et al. 2006). 



We conclude that the lack of the optical counterpart, the X-ray spectral properties and 
possible extended emission surrounding the source are all consistent with the proposed PWN 
identification. If the nebula was formed by the central pulsar associated with G27.8+0.6, its 
velocit}]^] of around 100-300 km s _1 would be consistent with measured velocities of other 
pulsars. In addition, its position close to the host SNR edge, suggests that perhaps the 
pulsar's velocity is at the transition from being subsonic in the SNR interior to supersonic 



outside the SNR. According to the current models (van der Swaluw 2005, and references 



8 this a factor of 2 below the flux listed in Table |3| due to using different model parameters 

9 The velocity is estimated from the remnant's age and the traveled distance of ~12 pc if the pulsar's 
birth place is at the SNR center. 



Fig. 3.— Broad-band (0.7-10 keV) EPIC image of the PWN candidate (Src.l in Tables [3 
and[§ see also Figure Q combining the PN, MOS1, and MOS2 data (left) and only MOS1 
and MOS2 data (middle). The images are smoothed with a Gaussian of FWHM 16". The 
pixel size is 4". The intensity scale is linear from to 1.42838 cts s -1 pixel -1 . The two arrows 
marked (1) and (2) show two possible directions of motion of the pulsar candidate (see text 
for detail). A VLA image of G27.8+0.6 at 1.4 GHz (right) is extracted from MAGPIS 
(White et al. 2005). The sources detected in the XMM-Newton observations are overlaid 
on the radio map. The position of the PWN candidate (Src. 1) is marked by the same two 
arrows shown in the EPIC image. 
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Fig. 4.— Broad band (0.7-10 keV), binned MOS1+MOS2 (left) and PN (middle) images 
showing the region around the detected PWN candidate, which consists of a bright point 
source and possible extended emission. After screening for the periods of high background, 
the exposure times were 4.9 ks for PN, 10 ks for MOS1, and 8.3 ks for MOS2. The NIR 
J-band (2MASS) image of the same region (right) shows no bright sources that could be 
identified with the X-ray emission (see text). The circle with a radius of 10" marks the 
position of the source. 
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Fig. 5.— Left: VLA image of G28.8+1.5 at 0.325 GHz extracted from MAGPIS ( [White 
et al. 2005) probably showing part of the SNR shell. Right: Color-coded (red: 0.5—1 
keV; green: 1—2 keV; blue: 2—4.5 keV), exposure-corrected and smoothed EPIC image of 
the same region showing soft diffuse emission across the SNR and around its center. The 
cataloged X-ray sources are overlaid on both images. 



therein), the pulsar's velocity becomes supersonic with respect to the ambient gas when the 
pulsar approaches the regions of cooler material near the edge of the shocked ejecta (for a 
SNR in the Sedov phase this occurs at the ~ 2/3 of the distance from the SNR's center). 
In this evolutionary stage, there are two PWN components: in addition to the X-ray bright 
PWN enveloping the moving pulsar, there is also a relic PWN (center-filled radio plerion) 
formed in the supersonic SNR expansion stage and still emitting synchrotron radiation at 
radio frequencies, with significantly longer lifetimes. Although our X-ray and radio data 
seem to be entirely consistent with this picture, further observations are needed to confirm 
that the source is extended and to study its X-ray morphology in detail. 

Finally, we note that we cannot entirely rule out the possibility that source #3 is a 
NS associated with G27.8+0.6. This source is located close to the center of the plerionic 



radio emission and its candidate optical/NIR counterparts are ~6"/2" away (see Section 2.2) 



However, even in this case, source #1 would remain a PWN candidate without a radio SNR 



association, which is a rather common occurrence (e.g., see Kargaltsev & Pavlov 2008) 
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3.2. Supernova remnant G28.8+1.5 

The bottom right panel of Figure [T] shows an amorphous region (~10' in diameter) of 
extended X-ray emission around the central part of SNR G28.8+1.5. Although faint, the 
diffuse emission has apparently a soft spectrum (suggesting a thermal plasma) as it can be 
seen in the exposure-correctecf^j color-coded EPIC image (Figure [5j right). The comparison 
with the radio map (Figure |5j left) shows that the X-ray emission comes from the region 
inside the bright rim of radio emission (i.e., perhaps part of the radio shell). In addition to 
shell-type, center-filled and composite (a combination of shell and center-filled morphology) 
SNRs, there is around 8% of SNRs with radio shells filled with thermal X-ray emission, 
called mixed-morphology SNRs ( Rho fc Petre||1998 ). Although the mechanism of the central 



thermal X-ray emission is not well understood, almost all of the mixed-morphology SNRs 



appear to be interacting with clouds (e.g., Koo et al. 2001). While the XMM-Newton 



data suggest that G28.8+1.5 might be a mixed-morphology SNR, high spatial resolution 
observations are needed to see whether this emission can be resolved into possible point 
sources and truly diffuse emission coming from the SNR, and measure its spectrum. 

There are two point X-ray sources close to the center of the diffuse emission, source #1, 
which we identified as a star (or, more likely, a pair of unresolved stars), and source #5, 
whose optical and infrared counterparts suggest an AGN classification. However, due to the 
broad PSF of the EPIC cameras, a Chandra observation would be needed to confirm some of 
these identifications. In particular, we note that several sources (most of which are identified 
as AGNs) have large offsets from their candidate optical counterparts (e.g., sources #6, 8, 
9, 17, 18, 21, 22; see Table [6]), which put their identifications into question. Among these 
sources, particular attention deserves source #6, which coincides with the brightest radio 
emission region that could be a radio PWN. 



4. Conclusions 

We analyzed XMM-Newton observations of two center-filled SNRs, G27.8+0.6 and 
G28.8+1.5, in a search for their central pulsars/PWNe and/or other types of NSs asso- 
ciated with these SNRs. We compiled the X-ray source catalogs, measured the fluxes and 
hardness ratios, and also searched for cataloged optical, infrared and radio counterparts, 
which we used to classify these objects. 

We found that most of the detected sources in both fields are foreground stars. The 



To produce the exposure-corrected images, we used the exposure maps corrected for vignetting. 
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proposed classification of these sources is based on hardness of their X-ray spectra and optical 
and infrared properties of their candidate counterparts in the USNO-B1 and 2MASS catalogs. 
The remaining sources are classified as AGNs (perhaps also including several galaxies and 
CVs), although the separations between the X-ray and optical/infrared positions for these 
objects are significantly larger. Hence, some of the proposed optical counterparts may be 
coincidental (e.g., some of these sources might be NSs). 

We identified one possibly extended source (Src. #1 in G27.8+0.6) as a PWN candidate 
associated with G27.8+0.6. This bright source does not have an optical counterpart and its 
morphology and X-ray spectrum seem to be consistent with the proposed PWN classification. 
The position of the PWN candidate with respect to its host SNR suggests a possible scenario: 
the alleged pulsar is moving away from the SNR center with a velocity of 200—300 km s _1 , 
and its position is near the edge of the shocked ejecta, approaching the cooler regions of 
the ambient material. At this stage the pulsar's velocity becomes supersonic, and a bow- 
shock nebula is being formed. A high-resolution timing observation is needed to search for 
pulsations from this source. 

Source #3 in G27.8+0.6 and source #6 in G28.8+1.5 could also be NSs associated with 
the SNRs, based on their locations within the regions of brightest radio emission (radio 
plerions) combined with relatively low probability of optical/NIR associations. 

We also detected faint extended emission around the central region of G28.8+1.5 (ap- 
parently inside its partial radio shell). If this emission is associated with G28.8+1.5, it would 
put this source in a rare class of mixed-morphology SNR. 

This work was partially supported by NASA grants NNG05GR01G and NNX09AC84G. 
The publication makes use of data products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center /California Institute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation, and the USNOFS Image and Catalogue 
Archive operated by the United States Naval Observatory, Flagstaff Station 
( |http: //www.nofs.navy.mil/data/| fchpix/) . 
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Table 1: The XMM-Newton observation log 





SNR 


Obs. ID. 


Date 


Coordinates 
RA/DEC (J2000) 


EPIC PN 


EPIC MOS1 


EPIC MOS2 


G27.8+0.6 
G28.8+1.5 


0301880401 
0301880901 


2005-10-04 
2005-10-04 


18:39:49.9 -04:24:00 
18:39:04.9 -02:43:44 


4.93 (8.53) 
14.03 (14.93) 


10.16 (10.16) 
16.38 (16.56) 


8.36 (10.16) 
16.56 (16.56) 



a exposure time in units of ks after background screening 
b total exposure time in units of ks 
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Table 2: Count rate to energy conversion factors (ECF; absorbed fTux=count-rate x ECF) 
for the medium filter of the EPIC instruments in the energy bands B1-B5, assuming a power- 
law model with the photon index of 1.7 and Galactic absorption of (1.5 — 2) x 10 22 cm -2 
measured in the direction of G27.8+0.6 and G28.8+1.5 ( Dickey fc Lockman||1990" ). 



Detector 


Bl 


B2 


B3 B4 


B5 






(io- 12 


erg cm -2 ct _1 


) 


PN 


1.56 


2.07 


3.36 8.76 


28.31 


MOS 


6.22 


6.44 


8.19 19.88 


94.95 
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